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ABSTRACT 

Binaries are typically excluded from direct imaging exoplanet surveys. However, the recent findings 
of Kepler and radial velocity programs show that planets can and do form in binary systems. Here, 
we suggest that visual binaries offer unique advantages for direct imaging. We show that Binary Dif¬ 
ferential Imaging (BDI), whereby two stars are imaged simultaneously at the same wavelength within 
the isoplanatic patch at high Strehl ratio, offers improved point spread function (PSE) subtraction 
that can result in increased sensitivity to planets close to each star. We demonstrate this by observing 
a young visual binary separated by 4" with MagAO/Clio-2 at 3.9 /im, where the Strehl ratio is high, 
the isoplanatic patch is large, and giant planets are bright. Comparing BDI to angular differential 
imaging (ADI), we find that BDPs Scr contrast is ^ 0.5 mags better than ADPs within ^ 1" for the 
particular binary we observed. Because planets typically reside close to their host stars, BDI is a 
promising technique for discovering exoplanets in stellar systems that are often ignored. BDI is also 
2-4 X more efficient than ADI and classical reference PSE subtraction, since planets can be detected 
around both the target and PSE reference simultaneously. We are currently exploiting this technique 
in a new MagAO survey for giant planets in 140 young nearby visual binaries. BDI on a space-based 
telescope would not be limited by isoplanatism effects and would therefore be an even more powerful 
tool for imaging and discovering planets. 

Subject headings: instrumentation: adaptive optics — techniques: high angular resolution — stars: 
individual (HD 37551) — binaries: visual — planetary systems 


1. INTRODUCTION 

A fundamental limit to directly imaging planets on so¬ 
lar system scales is the “speckle floor.” This is the region 
extending out to ^ 175 from the star that is dominated 
by speckles: time-varying aberrations in the point spread 
function (PSE). Because the locations and brightnesses 
of the speckles change slowly with time, it is difficult 
to fully subtract them, even with the help of adaptive 
optics (AO). Speckle noise thu s dominates over p hoton 
and sky noise close to the star (jRacine et al .11199^ . lim¬ 
iting our ability to directly detect faint planets at small 

separations. __ 

Angular differential imaging (ADL lMarois et a02QQ6[ ). 
in which the target star serves as its own PSE reference, 
has helped improve contrast close to the star, but only 
to a certain level and usually at the cost of attenuating 
the fluxes of the off-axis companions. ADI also requires a 
sufficient amount of sky rotation to occur throughout the 
observations: the less the sky rotates, the more a com¬ 
panion’s flux at a given separation is attenuated. This 
means that planets located in the speckle floor region will 
be more attenuated than those farther out, making their 
detection even more challenging. 
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Telescopes located at Las Campanas Observatory, Chile. 

^ Department of Terrestrial Magnetism, Carnegie Institute of 
Washington, 5241 Broad Branch Road, NW, Washington, DC 
20015, USA; email: trodigas@carnegiescience.edu 

^ Department of Physics & Astronomy, University of 
Rochester, Rochester, NY 14627, USA 

^ Steward Observatory, The University of Arizona, 933 N. 
Cherry Ave., Tucson, AZ 85721, USA 
^ Hubble Fellow 
^ NASA Sagan Fellow 


Simultaneous differenti al imaging (SDI; iMarois et al.l 
I2QQ3I: iRacine et al.l 1 1 9991) can in theory completely re¬ 
move speckles because the reference PSE is the target 
itself imaged simultaneously at different wavelengths, 
but it requires the faint companions to have unique and 
prominent spectral features over narrow wavelength re¬ 
gions compared to the stars-an assumption that is so far 
uncertain. Eurthermore, the different Strehl ratios of the 
two simultaneously-acquired images preclude true pho¬ 
ton noise limited PSE subtraction. Eor these reasons, 
it has been difficult for current and past direct imaging 
surveys to detect low-mass planets at small separations. 

At wider separations, a handful o f planets have been 
direc tly imaged ( e.g., HR 8799 bcde. IMarois et al.ll2QQ8L 
IMol: Beta Pic b , jransfe et al.ll2009ll2010l: HD 95086 
b. lRameau et alJl2fll,4 . However, with such a small sam- 
pie size, it is difficult to make statistical inferences on 
the general properties of giant planets, and even harder 
to constrain formation and evolution theories. To allevi¬ 
ate this problem, the number of imaged exoplanets must 
increase. 


We know from several imaging survey s that giant plan- 


ets are rarelv found at wide seoarations (iNielsen & Glose 

2010 

:lWahhai et al.ll2013l:lJanson et al.ll2012l:lBiller et al. 

2013 

: iMaire et al.l 20141). We also know from radial ve- 


locity (RV) and transit surveys that giant planets do re- 
side at sma l l separations (semimajor axis a < 5 AU; 
Butler et al.l I2QQ6I: I Gumming et al.l I2QQ8I: iBorucki et al.l 
2Q11[ ). Therefore future imaging surveys must improve 
detection sensitivity close to the star. 

Detecting planets at close separations-without build¬ 
ing new telescopes/instrumentation-calls for a different 
PSE subtraction technique. To avoid self-subtraction 
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(which plagues ADI), the target star cannot be used as its 
own PSF reference. To avoid atmospheric requirements 
in the companion (which plague SDI), the target and 
PSF reference must be imaged at the same wavelength. 
“Classical” PSF subtraction, wherein observations of a 
reference PSF are interleaved with observations of the 
target, is inefficient (two or more stars imaged to search 
for planets around one) and suffers from a time-varying 
PSF. Therefore the only remaining option is to observe 
two (or more) stars at the same time at the same wave¬ 
length. 

This is not a new idea. Speckle holography has 
been demonstrated as a viable technique to infer the 
presence of faint companions usin g the PSFs of mul¬ 
tiple simultaneously imaged stars (jLohmann fc Weigeltl 
rr^lWeigeVillma lK ores HHqoI), as long as they are 
within the isoplanatic p atch at the observed wavelength 
(|Christou fc Heg^ Il986f) . In space, the Hubble Space 
Telescope has been used to search for circums tellar dust 
in a binary by using each compo nenPs PSF (iLow et al.l 
Il999[ ). From the ground with AO. lMcElwainl (j2QQ7[ ) used 
Keck/OSIRIS to search for companions in several close 
visual binaries (no companions were found). 

We can improve on these previous studies in three 
ways. First, we can take advantage of modern AO sys¬ 
tems, which now offer high Strehl ratio imaging, espe¬ 
cially at long wavelengths. This improves our ability to 
detect faint planets close to their host stars. Second, 
we can image near 4 /rm. At this wavelength, gas giant 
planets are though t to be bright (jBurrows et al.l [2QQ1 : 
iBaraffe et al.l l2QQ3f ) , and the isoplanatic patch is large 
(^ 10-30"). A larger isoplanatic patch means binaries 
with larger projecte d separations will h ave nearly iden¬ 
tical PSFs. Indeed, iHeinze et al.l (j2010f ) imaged several 
binaries at 4-5 /im with the MMT deformable secondary 
AO system and showed that subtracting each PSF from 
the other could be more effective than ADL This is pos¬ 
sible because the target is no longer used as its own PSF 
reference, so the flux from the companion will never be 
attenuated by self-subtraction. Furthermore because the 
PSF is imaged at the same wavelength as the target, the 
companion is not required to have unique spectral fea¬ 
tures relative to the star (as with SDI). 

Finally, we can take advantage of advanced post- 
pro cessm£_aJ£oritlnn^liTe_pr^ component analysis 

(PCA, iSoummer et al.l[2Q12[ ). which has been shown to 
significantly improve PSF subtraction. This raises the 
possibility of completely removing the speckle floor from 
each star during the data reduction stage and thus reach¬ 
ing the photon noise limit. The photon noise limit is a 
critical threshold because, once reached, longer integra¬ 
tions result in increased sensitivity. This is not the case 
with ADI, where at small separations speckle noise domi¬ 
nates and adequate sky rotation, rather than integration, 
is preferred. 

It is now clear that planets can and do fo rm in bina¬ 
ries (|Buchhave et al.l 1201 ll: iNgo et al. ] |2Q15[ ). With the 
added benefit of potentially superior PSF subtraction, 
visual binaries should no longer be excluded from direct 
imaging planet surveys. In this paper, we show that Bi¬ 
nary Differential Imaging (BDI), in which two stars are 
imaged simultaneously, at high Strehl ratio, at the same 
wavelength, and within the isoplanatic patch, both im¬ 
proves detection sensitivity at small separations by ^ half 


a magnitude and is 2-4 x more observationally efficient 
than current ground- and space-based direct imaging sur¬ 
veys, respectively. In Section 2, we describe the BDI 
technique. In Section 3, we present proof of concept data 
obtained at the M agellan Clay telescope using Mag AO 
(iClose et al.l 120101) and the nea r-infrared (NIR) camera 
Clio-2 (|Sivanandam et al1l2QQ6f ). In Section 4 we com¬ 
pare the achieved contrast of ADI and BDI. In Section 

5 we discuss the technical and scientific benefits of BDI, 
on both ground- and space-based telescopes. In Section 

6 we summarize and conclude. 

2. BINARY DIFFERENTIAL IMAGING: EXPEGTED 
SIGNAL-TO-NOISE 

Two stars are simultaneously imaged on a detector at 
the same wavelength with AO. The separation of the 
stars is within the isoplanatic patch at that wavelength. 

_ _ _ M — M jQ 

Star A is brighter than Star B by / = 10 , where 

M is the magnitude of a given star at the chosen wave¬ 
length. Consider two pixels (Pixel A and Pixel B) each 
having the same coordinates (radius, position angle) rel¬ 
ative to each star. Assume that a hypothetical planet 
resides around Star A and some of its flux falls in Pixel 
A; Star B has no planets. We will calculate the total 
signal to noise (S/N) of the planet’s flux that results 
from subtracting one pixel from the other. In real high- 
contrast imaging PSF subtraction, many pixels comprise 
an image that is subtracted from another image, but the 
case is more simply explained using a single pixel. 

The total flux in Pixel A is given by 

Ta = Ia^Ip^ S, ( 1 ) 

where I a is the flux from Star A, Ip is the flux from 
the planet, and S is the flux from the sky. The sky 
flux is removed by subtracting blank sky {S 2 ) so that 
Ta Ia ^ Ip ^ S — 82 - We will assume that the sky 
flux is completely removed so that the total flux in Pixel 
A is just the flux from the star and planet. The same 
is true for Pixel B (without planet flux). In both cases, 
however, both read noise {R) and photon noise from the 
star (cT/) and sky {(Jsky) contribute to the uncertainty. 
We will assume that the photon noise from the planet 
{(Tip) is negligible. Therefore the total noise in Pixel A 
is given by: 


(JTa = ^ 


IOTa^ 


OTa^ 


OIa dS 




- 

= ^/Ia+ 2 S + 2 R^, (2) 


where the stellar photon noise = ^/Ta^ and we have 
assumed that as = as 2 = o-gky = ^/S and that the read 
noise is the same in the target and sky images. The total 
noise in Pixel B is similarly given by 


o'Tb — \/ Ip ‘^8 -|- 2R? 

= \/1 a/ f + 25 ' + 2R^ , ( 3 ) 

where we have made the substitution Ip = Ia/ f- 

Scaling and subtracting the Pixel B flux from Pixel A, 
we then have: 

F = TA-fTB^ (4) 
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which = Ip if the PSFs are identical. However, this 
is an ideal case and is unrealistic. Instead, there will 
probably be some slight differences in the PSFs due to 
non-isoplanatism. Therefore the remaining ffux will be 
F = /p + i, where i is this (positive or negative) extra 
ffux. We can define a new variable G = F — i = Ta — 
fTp — i and calculate the total noise: 


= ^JIa +2S + 2i?2 + p{lA/f + 2S + 2i?2) + af 
= y^(l + f)lA + (1 + P){‘^S + 2i?2) + (j?, (5) 

where we have assumed that the uncertainty in the scal¬ 
ing factor <7/ ~ 0. In reality, this uncertainty will be 
non-zero, but it is typically much smaller than the other 
sources of noise and is therefore ignored here. 

The total S/N of the planet in the BDI case is given 
by 

S/NpDi — — ! ^ =' 

\/(l + f)^A + (1 + /^)(25' + 2B?) + cr? 

(6) 

Here we have made the substitution Ip = c Ia, where c 
is the contrast of the planet relative to the host star at 
the location of the planet. 

It is thus evident that as / increases beyond the ideal 
/ = 1, S/Nbdi decreases. In the case of Pixel B minus 
Pixel A, we again have Eq. [6] except Ia ^ Ib and / ^ 
1//. In this case, S/Nbdi is maximized as / ^ oc. If the 
two PSFs are truly identical, then = 0 and the noise 
in BDI PSF subtraction is limited by the brightness ratio 
of the two imaged stars, with better PSF subtraction (in 
terms of noise) occurring when the brighter star PSF is 
subtracted from the fainter star PSF. If the PSFs are 
not identical, then the noise is limited by the amount of 
non-isoplanatism (cr^). 


2.1. Comparison to ADI 

In the case of ADI PSF subtraction. Pixel B and Pixel 
A both correspond to Star A, but Pixel B is recorded 
on the detector some time after Pixel A. We will again 
assume that Pixel A contains the planet ffux, and that 
Pixel B contains no planet ffux. However, to account 
for the self-subtraction that is inherent to ADI, we will 
assume that Ip a{r)Ip^ where a is the attenuation 
factor that depends on the distance from the star r. 

In theory, the brightness ratio of Pixel A and B (/) is 
unity, but the PSFs will never be identical. Therefore 
the noise is limited by the differences in the PSFs. To 
illustrate, we can derive cfadi in the case where we have 
the same star imaged at two different times resulting in 
two different total fiuxes, Ta (which contains the planet 
ffux) and Tp- The subtraction of the two pixels gives 
F = Ta — Tp = a{r)Ip -1- k, where k is the residual 
ffux that results from the non-identical PSF subtraction. 
As before we can define a new variable H = F — k = 
Ta — Tp — k and calculate the noise: 





S = 50, IA = too, fj/e = 30, and c = 10“^. The S/N ratio also 
depends on the attenuation factor a(r), which affects only ADI. 
The dashed line denotes the division between BDI or ADI being 
preferred. As a(r) decreases, BDI is preferred even at large ai 
values. When a(r) is close to I, BDI will only win out if cr^ is close 
to CTfc. This leads to the expectation that in real data, BDI should 
supersede ADI close to the star. 

= 'sJ^IIa + 2 (2S -f- 2i?^) + (j^. (7) 

The total S/N of the planet in the ADI case is given by 

S/Nadi = . ^ (8) 

^/2IA + 2(2S + 2B2)+a^ 

Eq. [7] is identical to the ideal (/ = 1) BDI noise (Eq. 
ED, except for the different noise term ak (compared to 
ai). This term (speckle noise) can in practice be esti¬ 
mated by calculating the standard deviation of pixels in 
an annulus at a given radius; it typically increases closer 
to the star and can be orders of magnitude larger than 
the other noise sources (jRacine et al.lll99^ . Thus ADI 
will yield lower noise than BDI if ak < cr^. However, 
S/Nadi is also limited by the attenuation factor a(r). 
Fig. [T] shows the theoretical S/N for BDI (ideal, / = 1) 
and ADI, for different values of <7/^, and a(r). It will 
be our goal in this paper to determine with real data how 
these parameters affect the achieved S/N of each imaging 
technique. 

2.2. Comparison to Classieal PSF Subtraetion 

We can also compare BDI PSF subtraction to classi¬ 
cal reference star PSF subtraction (sometimes referred 
to as “Reference Differential Imaging,” or RDI), which 
is typically employed for space-based telescopes. We will 
assume we are observing with a coronagraph on the the 
James Webb Space Telescope (JWST) NIRCam instru- 
ment. We will not compare t o PSF “roll” subtraction 
(|Schneider fc Silverstonel[2QQl) because JWST’s limited 
roll angle (^ 10°) will likely preclude exoplanet imaging 
using this technique, at least for close-in planets. 

The noise due to BDI PSF subtraction with NIRCam 
(no coronagraplfl) is again given by Eq. [5l Eor simplic¬ 
ity, we will assume the ideal case of two identical binary 

^ Coronagraphs (in the image plane) cannot be used on NIRCam 
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stars (/ = 1), no sky or read noise, and no noise due to 
the instrument (i.e., distortion), so that (Tbdi = \/27a. 
We will also assume = 0, since there are no isopla- 
natism effects in space. The signal in the BDI case is 
simply cl A, so that: 


S / NBDI ^space — 


C I A 


(9) 


The noise due to coronagraphic PSF subtraction is sim¬ 
ilar to the BDI case (again assuming an ideal / = 1), ex¬ 
cept for the extra noise due to PSF subtraction residuals, 
cT/c- Therefore we have: (Jcoron = + 9^^^^ where 

g is the fractional decrease in ffux that is caused by the 
coronagraph. The signal of the planet is tc/^, where t 
is the throughput of the coronagraph. The total S/N is 
then: 


S/N, 


t c Ia 


coron,space 






( 10 ) 


Comparing the S/N terms, we can see that 
S / Nbdi , space ^ S / NQQrpQjigpQ^QQ when C]^ > 
g~^^/2IA{t‘^ — g)‘ This relation can also be expressed as 
<^/c > — g)(Jij^^ where is the stellar photon 

noise. For g 10“^ and t = 0.81 (|Krist et al.ll2QQ7[) . the 
S/N from BDI PSF subtraction is larger if (Jk ^ lOa/^. 
It is not yet clear what (Jk will be for JWST/NIRCam, 
but this threshold is a useful diagnostic. Even if (Jk 
is small (a few x the stellar photon noise), BDI may 
still be preferred because it can reach better IWAs than 
coronagraphic imaging at a minor cost of slightly larger 
noise. Furthermore, because BDI is essentially photon 
noise limited, with enough integration it can always 
supersede classical PSF subtraction. 


3. TESTING BDI: OBSERVATIONS AND DATA REDUCTION 

To test the feasibility of BDI-in particular, the sim¬ 
ilarity of two PSFs separated by several arcseconds 
on the sky-we carried out proof of concept observa¬ 
tions of a binary with MagAO/Clio-2. We observed 
the young visual binary HD 37551, which consists of 
a G7 (primary) and K1 (jTorres et al.l I2QQ6I1 separate d 
by ^ 4"at a distance of 77T9 pc (jvan LeenwenI I2QQ7I ). 
The primary (henceforth Star A) has a V magnitude = 
9.650T0.017 and th e secondary (hencefor th Star B) has 
V = 10.538T0.017 (jLampens et al.l[200Th . HD 37551 is 
thought to be a member of the AB Dor moving group 
(jTorres et al.l I2QQ6I: iMalo et ^l2Q13[ ). However, the Li 
depletion pattern of the low-mass AB Dor members and 
their color-magnitude diagram posi tions appear to be in¬ 
distinguishable with the Pleiades (jLnhman et al.l [200^ 
iBarenfeld et al.ll2Q13f). whose Li-depletion bounda ry age 
is 130T20 Myr (jBarrado v Navascnes et aI]l2QQ4l ). We 
therefore adopt this age for HD 37551. 

We observed the binary on the night of UT 26 Novem¬ 
ber 2014 at the Magellan Clay telescope of the Las 
Campanas Observatory (LCO) in Chile. We obtained 
diffraction-limited AO images (300 modes) of the bi¬ 
nary in the narrowband 3.9 /im filter (AA = 0.23 /im, 
Ac = 3.95 /im; hereafter [3.9]) using the narrow camera 


(plate scale = 0701585 per pixe0). Star A served as the 
guide star for the AO system. The instrument rotator 
was turned off to enable ADI; therefore Star B rotated 
with the sky around Star A. Observing conditions were 
good, with the seeing ^ 0775. Every few minutes, we 
nodded the binary off the detector to acquire images of 
blank sky. Neither star’s PSE was saturated in any of 
the recorded images. 

After discarding frames of poor quality (e.g., open 
loop), we acquired 189 images of the binary for a total of 
31.5 minutes of integration. The parallactic angle (PA) 
changed by 19.37° throughout the observations. While 
obviously not a large amount of sky rotation, this is suf¬ 
ficient for directly comparing the detection sensitivities 
of ADI and BDI. 

All data reduction was performed using custom scripts 
in Mat lab. We divided the images by the number of 
coadds, corrected them for non-linearity, and divided 
them by the integration time to obtain units of detec¬ 
tor counts/s. To remove the sky background, for each 
image of the binary we constructed an optimal sky im¬ 
age using all the sky images, such that the noise in a 
small rectangular box far from the stars was minimized. 
Next we determined the sub-pixel location of each star by 
calculating the center of light inside a 0725 aperture cen¬ 
tered on the stars. We then registered and cropped the 
images of each star. Einally, we corrected for bad pixels. 
Eig. [2] shows the median PSEs for Star A and Star B, 
respectively. As expected, the PSEs (and speckles) are 
nearly identical, despite the different stellar brightnesses 
and 4" separation. 

Next, we performed ADI and BDI PSE subtraction. 
To simplify the comparison of the two techniques, we 
performed PSE subtraction on only Star A. Note that 
according to Section [21 the PSE subtraction noise for 
BDI will be higher in this case, so our BDI results should 
be viewed as conservative. 

3.1. ADI Reduction 

We reduced the Star A data using our custom 
ADI-hPCA pipeline. Because the PSE remained in ap¬ 
proximately the same location on the detector (aside 
from the off-chip sky nods), the PSE was very sta¬ 
ble. This meant that fewer PC A modes were needed to 
achieve excellent PSE subtraction. Using fewer modes 
also helped to minimize companion attenuation, which 
was important here given the small amount of sky ro¬ 
tation (< 20°). We settled on using 2 modes (out of a 
possible 189) to construct the PSE for each image. No 
rotation requirements were applied (i.e., all target im¬ 
ages were available for PSE construction), because do¬ 
ing s o does not necessarily increase recovered companion 
S/N (jMeshkat et al.ll2QM )FI After PSE subtraction, we 
rotated the images by their PAs to North-up, East-left 
and median-combined them to make the final image. No 
bright off-axis point-sources were detected in the final 
image. 


3.2. BDI Reduction 

The BDI+PCA data reduction is in practice quite sim¬ 
ilar to the ADI+PCA reduction, with one difference: in- 


with BDI because they change the shape of one PSF but not the 
other(s) in the FOV. 


® http://zero.as.arizona.edu/groups/clio2usermanual/ 
® See Section [52] for a justification of this choice. 
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Fig. 2.— The HD 37551 binary imaged at [3.9]. (a). Median-combination of all the Star A images, (b) The same, but for the Star B 
images. Despite the brightness difference between the two stars (measured to be 0.86 mags at [3.9]) and the 4" separation, the PSFs are 
nearly identical. However, some small differences are evident that make the BDI PSF subtraction noise (ai) non-zero. The sky background 
residuals are also different, most likely due to variations in the detector arising from Star B rotating around Star A during the observations. 


stead of the target itself serving as the PSF reference for 
PC A, the Star B images are used. As with ADI, the 
number of PC A modes was again a free parameter. The 
number of modes for BDI should not a priori be identical 
to the number of modes used in the ADI reduction, since 
the PSFs used by ADI and BDI are not necessarily iden¬ 
tical. Furthermore, because ADI is already hampered 
by companion self-subtraction, fewer modes will be pre¬ 
ferred for ADI to maximize recovered companion S/N. 
For BDI, we found that while increasing the number of 
modes did increase the signal-to-noise (S/N) of recovered 
artificial companions, the increase was marginal and at 
the expense of computation time. Therefore we settled 
on 10 (out of 189) modes, which was a good compromise 
between recovered S/N and efficiency. We rotated the 
PSF-subtracted images by their PAs to obtain North-up, 
East-left and median-combined them into a final image, 
which was additionally unsharp masked to remove resid¬ 
ual background flux. No bright off-axis companions were 
detected in the final BDI image, in agreement with ADI. 

4. ADI VS. BDI: COMPARING CONTRAST LIMITS 

To fully test the capabilities of BDI, we computed the 
overall contrast limits achieved by ADI and BDI. We ac¬ 
complished this using a Monte Carlo approa ch similar 
to the one employed in iRodigas et al.l ()2Q15f ) whereby 
a single artificial planet with a random contrast (Am) 
at a random separation (r) and position angle (0) is re¬ 
peatedly inserted into the raw data and recovered using 
both ADI and BDI. We inserted and recovered a total of 
10,000 artificial planets, which were flux-scaled replicas 
of the Star A PSF. The allowed parameter ranges were 
r G [0'.'15, 1'.'55], 0 G [0, 27r), Am G [5, 9]. Fig. [3] shows 
examples of artificial planets recovered by both BDI and 
ADI. 

We computed the S/Ns of the recovered a rtificial plan¬ 
ets usi ng the updated definition of S/N from iMawet et all 
(I20TI) . which is appropriate for speckle-dominated re¬ 
gions close to the star. Fig. |4a| shows the separation 
and contrast values of planets that were detected with 


S/N > 5 by ADI and BDI. BDI generally detects plan¬ 
ets closer to the star than ADI, being ^0.5 mags bet¬ 
ter in terms of contrast within ^ 1". Another way of 
showing this is using a binned contrast map, wherein we 
counted the number of successful detections (S/N > 5) 
by both ADI and BDI in O'.'l x 0.25 mag boxes. Fig. 
130 shows that BDI tends to detect a higher fraction of 
planets closer to the star than ADI. In terms of planet 
masses , using the COND mass-luminosity evolutionary 
tracks (|Baraffe et al.l[200^ for a 130 Myr old star, a half 
magnitude improvement translates into ^ 1 Jupiter mass 
improvement in sensitivity for BDI over ADL 

5. DISCUSSION 
5.1. Benefits of BDI 

We have shown with 3.9 /im MagAO/Clio-2 data that 
when two stars are imaged simultaneously at high Strehl 
ratio within the isoplanatic patch, BDI can be used to 
potentially probe smaller planet masses at closer sep¬ 
arations than in the ADI case. In addition, BDI of¬ 
fers two unique technical advantages. First, PA rota¬ 
tion is no longer a concern because each component is 
used as the other’s PSF reference. This is a major im¬ 
provement with regard to observing efficiency and sky 
coverage and offers the observer much more flexibility 
in planning an observational sequence over the course of 
a night (although non-isoplanatism may become a con¬ 
cern at high airmass). Second, because BDI allows us to 
search both stars for companions simultaneously, bina¬ 
ries can be studied 2x (4x) more efficiently than with 
ADI (classical PSF subtraction). 

From an astrophysical perspective, BDI offers the 
unique advantage of targeting stars (binaries) that are 
not usually preferred in large direct imaging surveys 0 

^*^ IThalmann et al.l(|20I4l) are conducting a direct imaging search 
for planets in spectroscopic binaries. However, the detectable plan¬ 
ets in that survey are ones that orbit both stellar components (P- 
type orbits). Such planets are presumably distinct from planets 
that orbit single stars. 
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Fig. 3.— Examples of inserted and recovered artificial companions, (a-b): ADI- and BDI-recovered Am = 6 planets at r = l''25 and 
0 = 95°. (c-d): The same, but at r = 0^5 and 6 = 240°. In both cases, the S/Ns of the BDI-recovered planets are higher than in the 
ADI case. It is also evident that the BDI images contain significant residuals very close to the star, unlike the ADI images. This indicates 
that the binary PSFs were not perfectly identical. Nonetheless, BDI still outperforms ADI at small separations because ADI significantly 
attenuates companion flux. 


The S-type planets (planets orbiting a single component 
of the binary) that BDI would be sensitive to should 
form and re main mostly unhindered by distant stellar 
companions (jBonavita fc DesideralDOOfl) . Furthermore, 
we know that ^ half of “hot Jupiter” systems have stellar 
companions, and these companions seem to ha ve no effect 
on th e primary star-planet orbital alignment (|Ngo et al.l 
|2Q15[ ) . It is thus important to search for and detect plan¬ 
ets in binaries and test theoretical pr edictions for the for¬ 
mation and stability of su ch planets (jHolman fc WiegertI 
119991: [jang-Condellll2Q15[) . From a more general perspec¬ 
tive, because ^ half of all s tars reside in multiple sys¬ 
tems (|Raghavan et al.l |2Q1Q[) , it is important to search 
for planets aro und the other half of s tars that are typi¬ 
cally excluded (|Thalmann et al.ll2Q14D . 

As a first attempt at probing planet formation and 
evolution in wide binaries, we are conducting a BDI 
survey for young exoplanets in visual binaries using 
MagAO/Clio-2. To minimize noise (e.g.. Section [2]), we 


have selected binaries that have comparable brightnesses 
in the NIR (Am < 2)E3 Their projected separations are 
between ^ 2-10" such that they lie within the isoplanatic 
patch at 4 /rm, their PSFs do not overlap, and they fit 
on the Clio-2 detector (size 9"xl5"). The primary com¬ 
ponent of each binary has V < 15 so that the AO system 
can lock. The binaries are young (< 500 Myr) and close 
to Earth (< 150 pc). Based on these constraints, the sur¬ 
vey sample consists of ^ 140 binaries (280 total stars) in 
young associations in the southern hemisphere. 

5 . 1 . 1 . Potential Space-based BDI Survey 

This is not a strict requirement. As described in Section [2] a 
larger brightness ratio means that the noise will be smaller when 
the brighter star is subtracted from the fainter one, but it will 
also be larger in the opposite case. This essentially reduces the 
efficiency of the survey because sensitivity is better for one star 
and worse for the other. 
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Fig. 4.— (a): 5(7 contrast for artificial planets recovered by ADI and BDI. The green points denote planets that were detected at S/N > 
5 by ADI, while the blue points (most of which are hidden by the green points) denote those that were detected by BDI. BDI generally 
detects planets closer to the star than ADI and is ~ 0.5 mags better in terms of contrast within ~ 1". (b): Binned contrast map showing 
the frequency of successful detections by BDI compared to ADI. Yellow squares contain more S/N > 5 detections by BDI than ADI, and 
vice versa for the blue squares. Green squares contain equal numbers of detections by both methods. BDI tends to perform better closer 
to the star than ADI. 


BDI is not limited to ground-based telescopes with AO. 
The JWST NIRCam instrument will offer space-based 
NIR imaging with a large 2.2' FOV. This opens up the 
possibility of imaging very wide binaries-or even single 
stars that happen to be close to other field stars-greatly 
expanding the potential survey size. Furthermore the 
similarity of two simultaneously-imaged PSFs should be 
limited by the (likely small) variations in the detector 
(like distortion) rather than the atmosphere/AO correc¬ 
tion. 

The number of possible survey stars would be re¬ 
stricted by the saturation of the NIRcam detector and 
the amount of PSF overlap. Using the JWST NIRcam 
PSF simulator (webbpslP^. in the absence of noise 98% 
of the PSF fiux is contained with ^ 5". Therefore stars 
with separations between ^ 5-60" could be observed. 
Assuming bright stars {V < 6) saturate the detector, 
there are conceivably ^ several hundred young visual bi¬ 
naries and young single stars (located close to other field 
stars) that could constitute a deep and powerful exo¬ 
planet imaging survey. 

5.2. Limitations 

Here we discuss possible limitations of our study and 
the BDI technique. First, we tested the ADI/BDI com¬ 
parison as a function of achieved sky rotation, since ADI 
works better with larger rotations. We divided the data 
in half (94 images in the first half, 95 in the second half). 
The PA changed by a total of ^ 10° in each half, and the 
total integration time for each half was ^15 minutes. We 
injected and recovered artificial planets (as in Section |4]) 
and compared the recovered S/Ns to the S/Ns achieved 
using the entire data set (full rotation and integration). 
We found that while the S/Ns of ADI-recovered planets 
increased with more rotation, the S/Ns of BDI-recovered 
planets increased by similar (or larger) amounts. This 

http://www.stsci.edu/jwst/software/webbpsf 


is because BDI is essentially integration-limited, rather 
than rotation-limited, so doubling the integration time 
invariably increases the contrast. Furthermore, increased 
sky rotation also benefits BDI b y azimuthally-smearing 
out residuals (jMarois et al.ll2QQ^ . 

Second, we tested the use of a rotation threshold in the 
ADI data reduction. This would be expected to reduce 
companion self-subtraction and thereby possibly increase 
companion S/N; however, this would also mean that cr/c, 
the noise from speckles in PSF subtraction, would be 
larger due to the increased dissimilarity between the tar¬ 
get and reference PSFs. Therefore it is not obvious that 
the overall S/N should increase. To test this, we inserted 
1000 artificial planets (as in Section |4]) and recovered 
them with BDI and ADL However, in the ADI case, for 
each image we required that the PSF be constructed from 
images that had rotated by at least 1 full-width half¬ 
maximum (FWHM) at the companion’s location. We 
found that the inclusion of a rotation threshold did not 
significantly increase the ADI-recovered companion S/N. 

Third, we tested the combination of ADI and BDI in 
the data reduction to potentially increase detection sen¬ 
sitivity. Using our Mag AO data, we tried running BDI 
first, then ADI, and vice versa (with ADI on each star). 
We also tested making all images (Star A and Star B) 
available for PSF reconstruction, as well as the same but 
with a rotation threshold on the Star A images. Overall, 
achieved contrast was not improved compared to BDI 
alone. However, in other cases (different brightness ra¬ 
tios and/or binary separations), an ADI+BDI combina¬ 
tion may prove more useful than either technique alone. 

A limitation of BDI (that ADI is not affected by) is 
the possibility of companion fiux attenuation by planets 
around the other star. For example, if both stars have 
planets at similar separations and position angles, then 
the fiux of each will be attenuated by the other. How¬ 
ever, even if two planets have identical brightnesses, the 
final fluxes will only be partially attenuated because each 
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planet’s flux is smeared out by PA rotation during the 
observations. Furthermore, assuming no preference for 
planet locations in an image, the chances of two planets 
being at the same locations around two stars are ^2% 
at 0''15 and much smaller farther out. 

Finally, it should be noted that coronagraphs, which 
have become increasingly utilized in exoplanet imaging 
efforts, cannot be used in BDL Specifically, coronagraphs 
in the image plane cannot be used. This limits the sensi¬ 
tivity of any BDI survey. However, it would be possible 
to use the apodizing phase plate (APP.lKenworthv et al.l 
I2QQ7I: iSnik et aT]l2Q12l: iQtten et~ffll2Q14l) . a pupil-plane 
“coronagraph” that offers improved contrast close to 
the star. Because the device sits in the instrument 
pupil plane, it affects all stars imaged on the detector. 
MagAO/Clio-2 already has several APPs installed, which 
can be used to further increase planet detection sensitiv¬ 
ity with BDL 


6. SUMMARY 


We have shown binary stars should be considered vi¬ 
able future targets for large direct imaging surveys be¬ 
cause of the potential for improved PSF subtraction close 
to the stars using BDI. This is important because planets 
are thought to be rare at wide separations. Based on the 
results of our proof of concept BDI data, we are conduct¬ 
ing a new direct imaging search for exoplanets in young 
visual binaries with MagAO/Cho-2. This survey should 
be able to detect fainter planets closer to their host stars 
and will tell us about exoplanet demographics in stellar 
populations that have never been probed before. BDI on 
a space-based telescope like JWST/NIRcam would not 
be limited by isoplanatism effects and would therefore 
be an even more powerful technique for imaging and dis¬ 
covering exoplanets. 

We thank the anonymous referee for helpful comments 
and suggestions. EEM acknowledges support from NSE 
grant AST-1313029. J.R.M was supported under con¬ 
tract with the California Institute of Technology (Cal¬ 
tech) funded by NASA through the Sagan Eellowship 
Program. 
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